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AMENDMENTS TO THE CLAIMS: 

This listing of claims replaces aU prior versions and listings of claims in the plication: 
Testing of Claims : 

1. (Currently amended) A computer implemented method for predicting the 
structure of a memhrane-bound protein having a pluraUty of [[a-]]helical regions, comprising: 

providing an amino acid sequence for the membrane-bound protein; 

uj ii ig the amino acid scqucnoo to identjfyjngeae two ormoie ranges of amino acids ii^ 
th^ ^ni^nn ar.irt sequence as transmembikie regions of the membrane-bound protein; 

constructing ^:.rh nf two or mote helices in a set of helices for the transmembrane 
regions^ 

aad optimizing a helix bundle configuration for the set of helices using a first molecular 

dynamics simulation; ' 

after o ptmiizint^ helix bundle configuration, c onstructing pr^gormQre a plurality of 

inter-hcHcal loops to generate a fiJl-atom model of the membrane-bound protein; 

optimizing the full-atom model using a second molecular dynamics simulation and 
outputting a predicted structure for the tefflsmembraneioimd protein based on the 

second optimization. 

2. (Withdrawn) A computational model of the structure of a transmembrane protein 
having a pluraUty of a-helical regions, the computational model comprising: 

a computer-readable memory storing data describing an optimized predicted three- 
dimensional structure for Hie transmembiane protein, the optimized predicted structure being 
generated according to flie metfiod of claim I- 

3. (Currently amended) The method of claim 1, wherein: 



R(aiDAT1OI13B0(l45:18:1OPM[EasteniDaiflightTiFnel"SVR:^^^^^ 



10/13/2004 14:21 FAX 6508395071 FISH & RICHARDSON "2)004 

AppUcaat : NaganxjaiiVaidehi.ctal. Attorney's Docket No.: 066IJK0600I /CIT319I 

Serial No. : 09/816,755 
Filed : March 23, 2001 
Page : 3 of 13 



constructing "^f^^ nf two or more helices in t he set of helices for the transmembrane 
regions includes constructing ^ach of two or more fr^e^ canonical heUces corresponding to the 
transmembrane regions, calculating a minimum-energy configuration for each of the canonical 
heUces, and optimizing each of the canonical helices ^ oflscmbling a holbc bundle including co eh 
of the pot of holicoo. ond cnlculntinfl n m ini mum onorgy configuration for the helix bundle in a 
lipid bilayer. 

4. (Withdrawn) A computational method for modeling the structure of a 
transmembrane protein having a plurality of a-helical regions, the method comprising: 

providing amino acid sequence information and sequence aUgoment information for a 
transmembrane protein having a plurality of a-helical regions; 

using the amino acid sequence information and the sequence alignment information to 
predict a set of transmembrane segments of the transmembrane protein; 

constructing canonical helices for the predicted transmembrane segments and optimizing 
the canonical helices using a first molecular dynamics simulation; 

combining the optimized heUces based on the sequence aligranent information to form a 
helix bundle, and assembling the helbc bundle with a lipid bUayer to form a system helix bundle; 

optimizing the structure of the system helix bundle using a second molecular dynamics 
simulation; 

adding inter-heUcal loops to the system helix bundle to form a full atom model; 
optimizing the fiill atom model using a third molecular dynamics simulation; and 
outputting a predicted structure for the transmembrane protein based on the third 
optimization. 

5. (Withdravra) The method of claim 4, wherein: 
transmembrane protein is a G-protein coupled receptor. 

6. CWithdrawn) The method of claim 4, wherein: 
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an energy minimum is calculated for each of the canonical helices before forming the 
helix bundle. 

7. (Withdrawn) The method of claim 4, farther comprising: 

deteraiining the periodicity of hydrophobic residues identified in the amino acid sequence 
information; and 

identifying a plurality of lipid-accessible residues based at least in part on the identified 
periodicity. 

8. (Withdrawn) The method of claim 4, wherein: ; . 
combining the helices based on the sequence information to form a helix bundle includes. 

orienting the helbc axes according to the 7.5 A electron density map for rhodopsin, 

9. (Withdrawn) The method of claim 7, wherein: 

combining ibe heUces based on the sequence information to foim a helix bundle includes 
orienting the identified Upid-accessible residues to face the outside of the heUx bundle. 
10- (Withdrawn) Themethodof claim 4, wherein: 

the first molecular dynamics simulation is a torsional molecular dynamics simulation. 

11. (Withdrawn) TTie method of claim 4, wherein: 

the second molecular dynamics simulation is a rigid body molecular dynamics 
simulation. 

12. (Withdrawn) The method of claim 1 1 , wherein: 

the first molecular dynamics simulation is a Newton-Euler Inverse Mass Operator 
dynamics simulation. 

13. (Witiidrawn) The method of claim 4, wherein: 
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the third molecular dynamics simulation is a mixed mode molecular dynamics 
simulation. 

14. (\irithdrawn) Themethodof claim 4, wherein: 

at least the third molecular dynamics simulation includes a solvent approximation. 

1 5. (WiHidrawn) The method of claim 14, wherein: 
the solvent approHmation is a continuum solvation model. 

16. (Withdrawn) The method of claim 15, wherein: 

the. solvent approximation includes the Surfece Generalized Bom model. . 

17. (Withdrawn) The method of claim 15, wherein: 

the solvent approximation includes the Poisson-Boltanann description. 

18. (Withdrawn) The method ofclaim 14, wherein: 

the solvent approximation is an empirical approximation comprising estimating solvation 
j&ee energy as a fimction of solvent accessible protein surface area. 

19. (Withdrawn) The method of claim 4, vrfierein: 

the predirted structure is generated by performing the third molecular dynamics 
simulation for a time in the range from about lOOps to about 1 ns. 

20. (Withdrawn) A computational model of the structure of a transmembrane protein 
having a pluiaUty of a-heKcal regions, the computational model comprising: 

a computer-readable data storage medium storing data describing an optimized predicted 
three-dimensional stmcture for the transmembrane protein, the optimized predicted structure 
being generated according to the method of claim 4. 
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21. (WWidrawii) Thecomputationalmodelof claim 20, wherein: 
the transmerobrane protein is olfectory receptor S6. 

22. CWithdiawn) The computational model of claim 20, wiiereia: 
the transmembrane protein is olfectory receptor SI 8- 

23. (Withdrawn) Thecomputationalmodelof claim 20, wherein: 
the transmembrane protein is olfactory receptor S 1 9. 

24. (Withdrawn) Thecomputationalmodelof claim 20, wherein: 
the transmembrane protein is olfactory receptor S25. 

25. (Withdrawn) The computational model of claim 20, wherein: 
the transmembrane protein is olfectory receptor S46. 



26. (Withdrawn) The computational model ofclaim 20, wherein: 
the transmembrane proteui is olfactory receptor S50. 

27. (Withdrawn) A computer program product on a computer-readable medium for 
predicting the structi>re of a membrane-bound protein having a plurality of a-heHcal regions, the 
computer program product comprising instructions operable to cause a programmable processor 
to: 

provide an amino acid sequence for the membrane-bound protein; 

use the amino acid sequence to identify one or more transmembrane regions of the 

membrane-bound protein; 

construct a set of helices for the transmembrane regions and optimize a helix bundle 
configuration for the set of helices using a first molecular dynamics simulation; 
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construct a plurality of inter-heUcal loops to generate a fiill-atoiu model of the 

membrane-bound protein; 

optimize the full-atom model using a second molecular dynamics simulation; 

output a predicted structure for the transmembrane protein based on the second 
optimization. 

28. (Withdravm) A computer program product on a computer-readable medium for 
predicting the structure of a G-protein coupled receptor having a pluraUty of a-helical regions, 
the computer program product comprismg instructions operable to cause a programmable 

processor to: ' ' ' 

provide amino acid sequence information and sequence- aUgnment information for a G- 

protein coupled receptor; . 

use the amino acid sequence information and the sequence alignment information to 
predict a set of transmembrane segments of the G-piotein coupled receptor, 

construct canonical heKces for the predicted transmembrane segments and optimize the 
canonical helices using a first molecular dynamics simulation; 

combine the optimized helices based on the sequence alignment information to form a 
helix bundle and assemble the helix bundle with a Upid bilayer to form a system helix bundle; 

optimize the structure of the system helix bundle using a second molecular dynamics 
simulation; 

add inter-helical loops to the system helix bundle to form a foil atom model; 
optimize the fiiU atom model using a third molecular dynamics simulation; and 
output a predicted structure for the G-protein coupled receptor based on the second 
optimization. 

29. (Withdrawn) A computational model of the structure of a G-protein coupled 
receptor, the computational model comprising; 
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a computer-readable data storage medium storing data describing a three-dimensional 
structure for olfactory receptor S6, the data describing the predicted structure including structure 
coordinates for the amino acids of olfactory receptor S6 having a root mean square deviation 
from Ae structure coordinates of the backbone atoms of the amino acids as set out in Table 2 of 
less t^fjfi or equal to aboiit 2,0 angstroms. 

30. (Withdrawn) A computational model of the structure of a G-protein coupled 
leceptorj the computational model comprising: 

a computer-readable data storage medium storing data describing a three-dimensional 
structure for olfactory receptor SI 8, the data describing the predicted structure including 
structure coordinates for the amino acids of olfectory receptor S18 having a root mean square 
deviation ftom the structure coordinates of the backbone atoms of the amino acids as set out in 
Table 3 of less than or equal to about 2.0 angstroms. 

31. (Withdrawn) A computational model of the structure of a G-protem coupled 
receptor, the computational model comprising: 

a conrputer-readable data storage medium storing data describing a three-dimensional 
structure for olfactory receptor SI 9, the data describing the predicted structure including 
structure coordinates for the amino acids of olfactory receptor S19 having a root mean square 
deviation from the structure coordinates of the backbone atoms of the amino acids as set out in 
Table 4 of less than or equal to about 2.0 angstroms. 

32. (Withdrawn) A computational model of the structure of a G-protein coupled 
receptor, &e computational model comprising: 

a computer-readable data storage medium storing data describing a three-dimensional 
stractuxe for olfeotory receptor S25, the data describing the predicted structure including 
stracture coordinates for the amino adds of olfactory receptor S25 having a root mean square 
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deviation from the structure coordinates of the backbone atoms of the amino acids as set out in 
Table 5 of less than or equal to about 2.0 aixgstroms. 

33. (Withdiavoi) A computational model of the structure of a G-protein coupled 
receptor, the computationai model comprising: 

a computer-readable data storage medium storing data describing a three-dimensional 
structure fox olfectoxy rece^nor S46. the data describing the predicted structure including 
structure coordinates for the amino acids of olfactory receptor S46 having a root mean square 
deviation from the structure coordinates of the backbone atoms of the amino acids as set out m 
Table 6 of less than or equal to about 2.0 angstroms; 

34. (Withdrawn) A computational model of the structure of a G-protein coupled 

receptor, the computational model comprisii^: 

a computer-readable data stomge medium storing data describing a three^imensional 
structure for olfactory receptor S50, the data describing the predicted structure including 
structure coordinates for the amino acids of olfectory receptor S50 having a root mean square 
deviation from the structure coordinates of the backbone atoms of the amino acids as set out m 
Table 7 of less than or equal to about 2.0 angstroms. 

35. (New) The method of claim 3, wherein: 

optimizing a helix bundle configuration includes calculating ^ minimum-energy 
configuration for the helix bundle in a Upid bilayer. 

36. (New) The method of claim 1, wherein: 

the membrane-bound protein is a G-protein coupled receptor. 

37. (New) The method of claim 1, wherein: 
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identifying two or more ranges of amino acids in the amino acid sequence as 
transmembrane regions includes aligning the amino acid sequence with an experimental or 
theoretical helical template. 

38, (New) The method of claim 1 , wherein: 

identifying two or more ranges of amino acids in the amino acid sequence as 
transmembrane regions includes determining the periodi^^^ 

amino acid sequence; and 

optimizmg a helix bundle configuration includes identifying a plurality of lipid-accessible 

residues based at least in part on the detennined periodicity. 

39. (New) The method of claim 1; wherein: 

constructing each of two or more helices in a set of helices for the transmembrane regions 
includes optimizing each of the two or more heUces in the set of heUces using a torsional 
molecular dynamics method. 



40. (New) The method ofclaim 39, wherein: 

the torsional molecular dynamics method uses the Newton-Euler Inverse Mass Operator. 
41 » (New) The method of claim 1, wherein: 

constructing each of two or more heUces in a set of helices for the transmembrane regions 
includes determining 3-D coordinates tiiat define tiie structure of each helix in the set of heUces. 

42. (New) The method of claim 1 , wherein: 

optimizing a helix bundle configuration includes deteiminmg a rotation and tilt of each 
helix in the set of helices. 
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43 . (New) The method of claiin 1 , A^vherein: 

optimizing a helix bundle configuration includes orienting the helix axes according to the 
7.5 A election density map for rhodopsin. 

44. (New) The method of claim 3 8, wherein: 

optimizing a helix bundle configuration includes orienting Ihe identified Upid-accessible 
residues to fece the outside of the helix bundle. 

45. G^ew) The method of claim 1 , wherein: 

the first molecular dynamics simulation is arigid body molecular dynamics simulation. 

46. OSTew) The method of claim 1, wherein: 

optimizing a helbc bundle configuration for the set of heUces includes modeling the effect 
of the environment of the membrane-bound protein. 

47. ^ew) The method of claim 45, wherein: 

the fiist molecular dynamics simulation uses the DREIDING force field, charges 
simulating the membrane, and charges for the transmembrane protein. 

48. (New) The method of claim 1 . wherein: 

the second molecular dynamics simulation is a mixed mode molecular dynamics 
simulation. 

49. (New) The method of claim 48, wherein: 

the second molecular dynamics simulation uses a torsional molecular dynamics method 
to model the heUces and inter-heUcal loops and a rigid body molecular dynamics method to 
model the membrane of the transmembrane protein. 
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50. (New) The method of claim 1 , wherein: 

the second molecular dynamics simulation includes dynamic optimization of the structure 
using cell multipole methods or fast torsional dynamic methods. 

5 1 . (New) The method of claim 1 , wherein: 

at least the second molecular dynamics simulation includes a solvent approximation. 

52. (New) The method ofclaim 51, wherein: 

the solvent approximation is a continuum solvation model. 

53. (New) The method of claim 52, wherein: 

the solvent approximation includes the Surface Generalized Bom model or the Poisson-. 
Boltzmann description- 

54. (New) The method of claim 53, wherein: 

the solvent approximation is an empirical approximation comprising estimating solvation 
free energy as a function of solvent accessible protein surface area. 

55. (New) The method of claim 1, wherein: 

the predicted structure is generated by performing the second molecular dynamics 
simulation for a time in the range from about lOOps to about 1 ns. 



56. (New) The method of claim 1, wherein: 
the set of helices includes four or more helices. 

57. (New) The method of claim 1, wherein: 
the set of helices includes seven or more helices. 
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